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ABSTRACT 

Azimuthally asymmetric dust distributions observed with the Atacama Large Millime- 
ter/submillimeter Array (ALMA) in transition disks have been interpreted as dust traps. We 
present Very Large Array (VLA) Ka band (34 GHz or 0.9 cm) and ALMA Cycle 2 Band 9 (680 
GHz or 0.45 mm) observations at a 0"2 resolution of the Oph IRS 48 disk, which suggest that 
larger particles could be more azimuthally concentrated than smaller dust grains, assuming an 
axisymmetric temperature field or optically thin 680 GHz emission. Fitting an intensity model 
to both data demonstrates that the azimuthal extent of the millimeter emission is 2.3 ±0.9 times 
as wide as the centimeter emission, marginally consistent with the particle trapping mechanism 
under the above assumptions. The 34 GHz continuum image also reveals evidence for ionized 
gas emission from the star. Both the morphology and the spectral index variations are con¬ 
sistent with an increase of large particles in the center of the trap, but uncertainties remain 
due to the continuum optical depth at 680 GHz. Particle trapping has been proposed in planet 
formation models to allow dust particles to grow beyond millimeter sizes in the outer regions 
of protoplanetary disks. The new observations in the Oph IRS 48 disk provide support for the 
dust trapping mechanism for centimeter-sized grains, although additional data are required for 
definitive confirmation. 


Subject headings: instabilities — protoplanetary disks — planets and satellites: formation 
interactions 


- planet-disk 


1. Introduction 

Studies of transitional disks - protoplanetary 
disks with inner holes in their dust distribu¬ 
tion - are revolutionizing our un derst anding of 
plane t formation (see the review bv lEspaillat et al~ 
I2014I1 . High spatial resolution Atacama Large 
Millimeter/submillimeter Array (ALMA) obser¬ 
vations reveal not only the dust cavities, but 


also azimuthal asymmetries (Iva n der Mare et_ah 
2013 : Gasassus et al.l 120131 : Fukagawa et, al. 20131: 


Perez < 


et al.ll2014lh and gas still present inside the 


dust c avities (van der Marel et al. 20 IT ICasassus et al 


2013: Bruderer et al 


2014: Zhang et al.l 201 if 


van der Marel et alil2015f Perez et al.120151 ). Dust 


trapping has been s uggested as a solution for the 
radia l drift problem ( Whinnlell 1972 : Weidenschillina 


WENI fellow 


19771 : Brauer et ah 120081 ) which prevents dust 
grains from growing beyond millimeter-sizes in 
the outer disk. Trapping facilitates the crucial 
steps in dust growth toward the formation of plan- 
etesimals (and therefore planets). Because dust 
continuum emission is dominated by grains with 
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sizes up to t hree times the observing wavelength 
( Drainel 20061) , different continuum wavelengths 


probe different particle sizes. Since particle trap¬ 
ping depends on particle size, disk turbulence, and 
the pressure gradient profile, multi-wavelength ob¬ 
servations are required to confirm the trapping 
scenario and constrain disk param eters such as 
viscosity (e.g., Birnstiel et al. 20131 and references 


therein). 

The Oph IRS 48 transition disk exhibits a 
highly asymmetric structure in the 0.45 mm dust 
continuum. The continuum emission is generated 
by millimeter-sized grains gathered in a peanut¬ 
shaped structure, spanning less than a third of 
the disk ring azimuth; the peak of this structure 
is > 100 x higher than the opposite side of the 


is > iUUx_mgnei_tnan tne opposite side ot tne 
ring ( van der Marel et al • m 3. However, the gas 


traced by 12 CO <7=6-5 line emission has an ax- 
isymme tric disk distribution down to 20 AU in 


radius ([Bruderer e t al.1 |2014j), confirmed by CO 


isotopologue observations of the same disk (van 
der Marel et al. submitted). Likewise, the ther¬ 
mal mid-infrared and near-infrared scattered light 
emission, tracing small micrometer-sized grains, 
suggest a ring-like structure although variations 
along t he ring cannot be seen due to high op tical 
depth ( Geers et al.l 120071 : iFollette et al.l l201ai . A 


separation between large and small grains/gas in- 


trap 

Bare:e & Sommeria 

1995 

Klahr & Henning 

1997; 

Brauer et al. 200(1; 

Pinilla et al. 20121). Az- 

imuthal traps may result 
be due to instabilities si 
or baroclinic instability 

from vortices, which can 
ich as the Rossby wave 

(e.g., 

Lvra & Lin 201s 

Raettig et al. 2013; Fung et al. 2014[ Flock et al. 


20151). A vortex locally increases the gas pressure 
over a limited radial and azimuthal extent. In the 
outer disk, millimeter-sized particles will then drift 
toward this pressure maximum a nd get trapped in 


the azimuthal dire ction (e.g., lAtaiee et al.l 12013 : 


Zhu fe Stone 201 1;, which can explain the ob¬ 
served features in IRS 48. Inside dust traps, the 
dust particles continue to grow to the maximum 
grai n size permitted by the fragmentation bar¬ 
rier ( Birnstiel et al. 2010h , which in t he outer disk 
corre sponds to centimeter-sizes ('e.g.. lPinilla et af 


2012). Further growth to planetesimal sizes can 
occur by streaming instabilities (IJohansen et al 


2009h or b y tak ing mas s tran sfer effects into ac¬ 
count ( Windmark et alJ 2012]). 


The trapping depends on the particle size: 
small dust particles are strongly coupled to the 
gas and thus follow the gas distribution, while 
larger particles (traced at longer wavelengths) are 
less coupled and therefore feel the drag force to¬ 
ward the pressure maximum, resulting in a more 
spatially concentrated distribution of these par¬ 
ticles, both radially and azimuthally. Both the 
trapping concentration and the maximum particle 
size in the trap are depende nt on the turbulenc e 
and the gas surface density ( Pinilla et al. 20151) . 
In this paper, we present spatially resolved ob¬ 
servations of Oph IRS 48 at both millimeter and 
centimeter wavelengths, taken with ALMA and 
the Karl G. Jansky Very Large Array (VLA). The 
distributions of dust emission at 450 ^m and 0.9 
cm are compared and we aim to provide confirma¬ 
tion of the dust trap scenario in a transitional disk. 
Section 2 describes the observation setup and cal¬ 
ibration process, Section 3 presents the images 
and the derived intensity profile and Section 4 
discusses the interpretation of the different mor¬ 
phologies and possible implications for the dust 
trapping scenario. 


2. Observations 

Observations of Oph IRS 48 at 34 GHz (9 mm) 
were obtained using the VLA in the CnB and B 
configurations in January-February 2015 as part 
of program 14B-115, with baselines ranging from 
75 to 8700 m. The spectral windows were con¬ 
figured for a maximum possible continuum band¬ 
width of 8 GHz centered at 34 GHz (Ka band) 
in dual polarization, using 64 spectral windows of 
128 MHz each, with 3-bit sampling. Due to the 
low declination, the source was observed in three 
scheduling blocks of 2.75 hours with 1.25 hours 
on source in each block. The bandpass was cali¬ 
brated using J1517-2422 (in the third block J1924- 
2914 was used instead), the absolute flux was cal¬ 
ibrated using 3C286 and J1625-2527 was used as 
gain calibrator, periodically observed every three 
minutes. The pointing was checked on the gain 
calibrator in X-band every 30 minutes. The cali¬ 
brated data were concatenated and imaged using 
Briggs weighting with a robust parameter of 0.5. 
Deconvolution using CASA imfit reveals that the 
source is marginally resolved. The flux calibration 
uncertainty is 10%. 
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Table 1: Image properties 


Telescope 

Frequency 

(GHz) 

Bandwidth 

(GHz) 

Beam size 

Beam PA 

Flux 

(mJy) 

Peak 

(mJy beam -1 ) 

rms 

(mJy beam -1 ) 

ALMA 

680 

4.7 

0"19x0"15 

79° 

1000 

190 

0.6 

VLA 

34 

8.0 

0"46x0"26 

21° 

251xl0 -3 

138xl0 -3 

3.5xl0 -3 


F (mJy) 


F (/.*Jy) 


o 

<u 

Q 

< 


1.0 0.5 0.0 - 0.5 - 1.0 0.5 0.0 - 0.5 - 1.0 0.5 0.0 - 0.5 - 1.0 0.5 0.0 - 0.5 - 1.0 0.5 0.0 - 0.5 - 1.0 

ARA (") ARA (”) ARA (’’) ARA (”) ARA (”) 

Fig. 1.- ALMA and VLA observations of dust emission at 680 and 34 GHz of the Oph IRS 48 disk, (a) 
The 680 GHz image (color scale); (b) the 34 GHz (color scale); (c) the overlay of the 34 GHz contours in 
red (taken at 3,6,...,39er, cr=3.5 /iJy) on the 680 GHz contours in black (taken at 3,12, 24,...,324cr, with 
<r=0.6mJy); (d) the overlay of the observed 34 GHz image (red) on the 680 GHz image imaged with the 34 
GHz beam (black) in contours taken at 10%,20%,...,90% of the peak; and (e) the overlay of the 680 GHz 
image as if observed by the VLA (sampled on the VLA visibilities) with the 34 GHz observations taken at 
3,6,9,12,...er, with a = 3.5/tJy. Ellipses indicate the beam size listed in Table [Q 



Observations at 680 GHz (440 /mi) were ob¬ 
tained using ALMA in Cycle 2 in July 2014 in the 


C34- 3 configuration with Band 9 (jBarvshev et al 


2015h as part of program 2013.1.00100.S. The ob¬ 
servations were taken in four spectral windows of 
1920 channels: three windows have a bandwidth 
of 937.5 MHz, centered on 661, 659 and 675 GHz, 
the fourth spectral window was centered on 672 
GHz with a bandwidth of 1875 MHz. The total 
continuum bandwidth was ~4.7 GHz. The flux 
was calibrated using J1517-243, the bandpass with 
J1427-4206 and the phase with J1626-2951. The 
total on-source integration time was 52 minutes. 
The data were self-calibrated and imaged using 
Briggs weighting with a robust parameter of 0.5. 
The flux calibration uncertainty is 20%. 

Table □ lists the properties of the images. The 
astrometric accuracy is set by the calibrators, and 
is typically ^30 mas for these two data sets, which 
is much smaller than the beam size. 


3. Results 

Figure [T| shows images of the dust continuum 
emission at 680 and 34 GHz. The 680 GHz con¬ 
tinuum has a simila r asy mmetry as observed in 
van der Marel et al. ( 2013 ) in the Band 9 Cycle 0 
data, though these new data have improved angu¬ 
lar resolution. The 34 GHz continuum emission 
peaks at the same location, but is much more az- 
imuthally concentrated. This is not a sensitivity 
effect: when the 680 GHz data are restored with 
the same cleaning beam as the 34 GHz data, its 
azimuthal extent is still clearly wider than in the 
measured 34 GHz data (see Figure QJI). The ra¬ 
dial width cannot be compared due to the verti¬ 
cally elongated VLA beam shape, which is caused 
by the low declination of the source with respect 
to the VLA site, but given the disk geometry, this 
elongation does not affect constraints on the az¬ 
imuthal width, the key parameter of interest here. 
Spatial filtering can be ruled out as an expla¬ 
nation for the different azimuthal extents; these 
data recover all of the flux found on the short- 
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est baselines from previous VLA observations in 
the DnC configuration (beam size 3.3x 1.3”, flux = 
252 ±ll/iJy). Also, we have simulated the ALMA 
image as if observed by the VLA by sampling the 
ALMA image on the VLA visibilities (Figure QJ) 
to rule out horizontal spatial filtering. Besides, the 
largest angular scales recovered by the CnB obser¬ 
vation is 5”, which is much larger than the disk. 
Thus, the centimeter-sized dust grains (traced by 
the 34 GHz continuum) have a narrower azimuthal 
distribution than the millimeter-sized dust grains 
(traced by the 690 GHz continuum). The peak 
brightness temperatures for the 680 GHz and 34 
GHz are 31 and 1 K, respectively. 


In order to quantify the concentration of the 
dust grains, the morphology of the continuum im¬ 
age is fit to a two-dimensional fourth power Gaus¬ 


sian intensity profile I(r,4>) (jv an d er Marel et al 


201311 : 


I{r, <j>) = I c exp 


~(r ~ r c ) 4 

2 ri 


exp 


2 # 


(!) 

Previous studies describe azimuthal dust asymme¬ 
tries d ue to vortices as regular second powe r Gaus- 
sians ( Lyra &; Lin 20131: Perez et al. 2014 ). but it 
was found that our data are fitted much better 
with the A th power equation. The stellar posi¬ 
tion was set to 16 ?l 27 m 37 s .185 -24°30’35"39, taken 
from a Keplerian model fit to the CO isotopologue 
data from the same ALMA data set (van der Marel 
et al. submitted). The proper motion is negligi¬ 
ble compared to the beam size for the period in 
between these observation sets (~6 months), so 
they can be overlaid directly. For the fitting, we 
use x 2 minimization over a grid of </> c , </>w, U: and 
r w with steps of 1° and 1 AU for angle and ra¬ 
dius respectively. The fitting is performed in the 
uv-plane, with the model visibilities sampled onto 
the observed spatial frequencies. We find the best 
fit (see Figure [2]) for r c =61(± 2) AU, r w = 14(±2) 
AU, <^ c =100(±3)° with <f> w = 41(±4)° for the 680 
GHz and (j> w =18(±7)° for the 34 GHz continuum. 
The area at 680 GHz is smaller than that reported 


van der Marel et al.1 (l2013i) . due to better spa¬ 


tial resolution, although still within the derived 
error bars. The values and errors for r c , r w and 
(j> c are mainly constrained by the ALMA data, as 
the VLA data did not have sufficient spatial reso¬ 
lution to constrain the error bars. The errors are 
estimated by rescaling the reduced \' 2 to 1. 
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Fig. 2.— Best fits of the intensity profile given 
in Equation (1) to the dust continuum for the 680 
GHz (left) and 34 GHz (right) data. The top panel 
shows the data in colorscale (the yellow contour 
shows the 3er level), the middle two panels the 
unconvolved and convolved model, and the bot¬ 
tom panel the residual image, with the absolute 
values of the residual in the same color scale as 
the data image, and overlaid contours at 5er spac¬ 
ing (ct 6 80GHz=0.6 mjy beam -1 , <j 3 4GHz=3.5 fiJy 
beam -1 ). Dashed contours indicate negative in¬ 
tensity levels. 
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The important conclusion is that the azimuthal 
width of the centimeter emission is 2.3(±0.9) times 
narrower than the millimeter emission. The aspect 
ratio of the submillimeter emission is 3.1 (±0.6), 
and >1.4 for the centimeter emission. 

The residual of the 680 GHz image still shows 
significant emission outside the Gaussian fit (peak 
25er), especially at the tail in the east: 15% of 
the total absolute flux remains in the residual. 
The intensity equation described by Eqn. [T] is 
clearly not sufficient to describe the detailed struc¬ 
ture: additional (vertical) features may be present. 
The disk is known t o have a large scale height 
( Bruderer et al . m , which is not taken into ac¬ 
count in the intensity model, while the dust is pos¬ 
sibly optically thick at this wavelength and the 
vertical structure may be relevant. However, as 
still more than 85% of the structure is recovered 
by the Gaussian intensity equation and the SNR is 
very high, the description is sufficient to compare 
the azimuthal width at the two wavelengths. 

The residual of the 34 GHz image clearly shows 
a point source at the stellar position, with a peak 
flux of 36 fiJy beam -1 (lOcr). This point source 
can be either dust emission from an unresolved 
inner disk, or free-free/synchrotron emission from 
ionized gas close to the star. Since this emis¬ 
sion is not seen in the ALMA data, the spectral 
index of this emission is <1.3, so dust emission 
is unlikely. The origin of this emission can be 


determined usin 

g long 

er wavelength observations 

(Rodmann et al. 

2006 

Ubach et al. 2012). The 


total flux of the dust emission in the dust trap 
is thus only 216 /zJy at 34 GHz. 


4. Discussion and summary 

The difference in azimuthal width between the 
two continuum images suggests a segregation of 
particle sizes, where centimeter-sized grains are 
more azimuthally concentrated than millimeter¬ 
sized dust grains. However, optical depth effects 
can hide a narrow concentration of millimeter-dust 
grains equivalent to the 34 GHz morphology on 
top of the apparent millimeter-dust distribution. 
The optical depth r 680 Gi?z cannot be measured in¬ 
dependently from the temperature with the avail¬ 
able data. The measured brightness temperature 
at 680 GHz of 31 K (assuming the Planck equa¬ 
tion) is only a factor of two lower than the cal¬ 


culated physical dust temperature of 60 K (at 60 
AU r adius in a physical disk model. iBruderer et al.l 


2014) so the emission is likely not highly optically 


thick. In order to quantify the optical depth at 680 
GHz, we assume that the peak millimeter emission 
originates from an isothermal area flbeam (sr) of 
the beam size with Td U st=60 K. This is a conserva¬ 
tive limit since r likely decreases azimuthally out¬ 
ward. For the peak flux of 189 mJy, T68OGffz=0.5, 
so the emission is indeed marginally optically thick 
at 680 GHz. If the dust temperature is much 
lower than assumed here, the emission becomes 
fully optically thick and a narrow concentration of 
millimeter-dust grains could remain hidden. Lo¬ 
cally lowering the temperature is possible by shad¬ 
owing due to an incli ned inner disk , such as pro¬ 
posed for HD142527 ( Marino et al.l l2015l) to ex¬ 
plain the scattered light emission. For Oph IRS 48, 
there is no evidence for a local temperature drop, 
so we propose that the difference in emission origi¬ 
nates from a spatial segregation between particles, 
where the centimeter-sized grains are more con¬ 
centrated than the millimeter-sized grains. Spa¬ 
tially resolved continuum observations at interme¬ 
diate wavelengths are required to confirm the dust 
temperature and optical depth. 

Combining the total 34 GHz flux of 216±22 
/ 1 Jy with SMA and disk integrated ALMA obser¬ 
vations at 230, 345 and 680 GHz gives a spectral 
slope of a = 2.84 ± 0.06(1% ~ v a ). The fluxes 
are 50±7.5, 160±24 and 1000±200 mJy respec¬ 
tively, taken from Brown et al.l ( 20121 ) and this 
work, with the errors dominated by the flux cal¬ 
ibration uncertainty. Measuring a independently 
between 680 and 230 GHz (2.8±0.2) and between 
230 and 34 GHz (2.85±0.7) results in the same 
value within error bars, supporting at most mod¬ 
est optical depth at 680 GHz. 

Because both the VLA and the ALMA data 
are spatially resolved, the spectral index can be 
computed as a function of position (see Figure [3]). 
For this figure, the ALMA data were convolved 
with the same synthesized beam as the VLA ob¬ 
servations. The figure shows that a increases from 
2.6±0.08 in the center to 3.3±0.15 in the outer 
wing s, similar to the synth etic spectral index map 
from Birnstiel et al. ( 20131) of azimuthal trapping. 
The error in spectral index across the image is 
calculated based on the signal-to-noise ratio, as 
shown in the right panel of Figure [3] Therefore, 
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Spectral index a Error spectral index a 


2.0 2.3 2.7 3.0 3.3 3.7 4.0 0.00 0.04 0.08 0.12 0.16 0.20 



1.0 0.5 0.0 - 0.5 - 1.0 0.5 0.0 - 0.5 - 1.0 

Relative RA (") Relative RA (") 


Fig. 3.— Comparison of the ALMA and VLA data 
of Oph IRS 48, imaged with the same synthesized 
beam of 0"46x0"26 (lower left corner). The left 
plot shows the spatially resolved spectral index a, 
using 34 GHz and 680 GHz data. The right plot 
shows the error on a as function of position. The 
ALMA data are overlaid in yellow contours with 
10%,20%,...,90% of the peak flux. The spectral 
index decreases at the peak of the dust emission. 
Data points where the 34 GHz flux density is less 
than 2a (a = 3.5/iJy) are excluded. 


the decrease of the spectral index at the peak of 
the dust emission is significant. Due to the flux 
calibration uncertainty, which is 20% for ALMA 
Band 9 and 10% for the VLA Ka band observa¬ 
tions, the calibration error over the entire image 
on a a is ~ 0.07, but this does not vary with po¬ 
sition and therefore it is not included in the error 
plot. 


Analysis of the dust opacity index /3 (t„ oc v ^) 
is an important tool in constraining the dust par¬ 
ticle sizes in the disk. When the dust emission 
is optically thin (r <C 1) and in the Rayleigh- 
Jeans regime, /3 = a — 2, where /3 < 1 indicates 
dust growth t o lar ger than millimeter-sized par¬ 
ticles ( Drainel [20061) . For IRS 48 the 680 GHz is 
not optically thin and these assumptions cannot 
be made. Therefore, the dust opacity index thus 
has to be derived from the optical depth r„ itself. 


Using r„ oc , /3 can be calculated as a function 
of position by calculating r„ at both frequencies at 
different positions, assuming Tdust=60 K and the 
flux filling the beam. We find /3 ~ 0.7 ± 0.1 in the 
center and /3 ~ 1.3 ± 0.3 in the outer parts for a 
uniform temperature. As the continuum peak is 
likely more optically thick in the center, the emis¬ 


sion there traces higher vertical layers where the 
temperature is higher, implying an even lower /3: 
for T dust =200 K in the center, /3 = 0.6 ± 0.1 ac¬ 
cording to the same calculation. The azimuthal 
trend in (3 is thus consistent with increased dust 
growth (increase of a max ) or an increase of larger 
particles in the center (increase of a m ; n ). On the 
other hand, if the temperature is in fact as low 
as 31 K, /3 would increase in the center. In or¬ 
der to get a uniform /3 along the entire azimuthal 
shape (which implies no change in dust growth or 
size segregation), Td us t needs to be as low as 24 
K at the center, which is more than a factor of 
two lo wer than the derived m idplane dust temper¬ 
ature ( Bruderer et al.1 [2014 1. while the optically 
thick emission likely originates from higher ver¬ 
tical layers with even higher temperatures. Thus, 
within reasonable azimuthal variations of the tem¬ 
perature, the results hint at variation in /3, though 
some uncertainties remain. 

In summary, the centimeter emission provides 
further support for the dust trapping mechanism 
in the Oph IRS 48 disk: within the assumptions 
of the 680 GHz optical depth and the tempera¬ 
ture field, the centimeter-sized grains appear to be 
more concentrated in azimuth than the millimeter¬ 
sized grains as predicted in analytical dust mod¬ 


els of azimuthal pressur e maxima (IBirnstiel et al 
2013 : Lvra fc Lii] 20131 ). Inside the dust trap, 


grains may have reached centimeter sizes (and per¬ 
haps even larger), the first step toward planetesi- 
mal and planet formation. 
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